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Abstract 

We compute the semi-classical quantum amplitude to go from an initial 
spherically symmetric bosonic matter and gravitational field configuration 
to a final radiation configuration, corresponding to the relic Hawking ra- 
diation from a non-rotating, chargeless black hole which evaporates com- 
pletely. This is obtained via the classical action integral which is solely a 
boundary term. On discretising the classical action, the quantum ampli- 
tude can be expressed in terms of generalised coherent states of the har- 
monic oscillator. A squeezed-state representation is obtained by complex- 
ifying the proper time separation T at spatial infinity between the initial 
and final space-like hypersurfaces. Such a procedure is deemed necessary 
as the two-surface problem for Dirichlet boundary data and wave-like per- 
turbations is not well posed. We find that infinitesimal rotation into the 
lower complex T plane is equivalent to a highly-squeezed final state for 
the relic radiation, similar to the relic gravitational-wave background in 
cosmology. This final state is a pure state, and so the unpredictability 
associated with the final momentarily-naked singularity is avoided. The 
cosmological analogy is the tunnelling from an initial smooth Euclidean 
or timeless state to a classical universe. The high-squeezing limit cor- 
responds to a final state of the Hawking flux which is indistinguishable 
from a stochastic collection of standing waves. The phases conjugate to 
the field amphtudes are squeezed to discrete values. We also discuss the 
entropy of the final radiation in the high-squeezing limit. 
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1 Introduction 



Much has been achieved in the apphcation of Feynman's quantum-amphtude for- 
mahsm to black-hole evaporation. In a semi-classical approximation, Hartle and 
Hawking [1] pioneered this technique, proving, by analytic continuation in the com- 
plexified space-time, that the amplitude for a black hole to emit a particle was 
proportional to the amplitude for the black hole to absorb a particle. The pro- 
portionality factor is just a Boltzmann factor, whence one determines that a black 
hole of initial mass M/ can come into equilibrium with the exterior radiation only 
at the Hawking temperature Th = hc^ /SnGMjkB- Gibbons and Hawking [2] again 
employed semi-classical path-integral methods to re-derive Hawking radiance and 
the Bekenstein-Hawking entropy formula. 

In this paper, we determine the semi-classical wave functional of specified bound- 
ary data in the context of classical-like coherent states and purely quantum-mechanical 
squeezed states. The emphasis will be on the final radiation which remains after 
a black hole has evaporated completely, and its comparisons with the relic Cosmic 
Microwave Background Radiation (CMBR) induced by cosmological perturbations. 

It is well-known that particle creation by black holes has many similarities with 
cosmological particle creation, despite the lack of asymptotic fiatncss in the cosmo- 
logical scenario. Cosmological and black-hole particle creation both require a time 
dependence in the metric. The time-dependent background gravitational field (also 
known as the 'pump' field) is a result of the non-linearity of Einstein's equations, 
which is displayed through the coupling of the gravitational wave field, say, with 
the evolving background space-time. One is also aware of cosmological No-Hair 
Theorems, analogous to the black-hole well the presence of singularities, 

at some finite time in the past in the cosmological case, and hidden behind event 
horizons in the black-hole case. 

Rather than employ a time-dependent background metric, black-hole evaporation 
can be interpreted as a tunnelling process, whereby the strong gravitational field in 
the neighbourhood of the future event horizon promotes vacuum fiuctuations into 
real particles. Positive-energy particles are detected by observers at infinity, while 
negative-energy particles tunnel through the event horizon reducing the mass of 
the black hole. Eventually the black hole presumably disappears completely (we 
assume in this work that there are no black-hole relics) leaving only radiation and an 
almost fiat space-time. The singularity inside the black hole is momentarily naked 
prior to total disappearance [20], and is deemed to take with it the information 
about possible collaspe configurations which created the black hole. This is the 
information-loss paradox. 

In cosmology, near the start of infiation, the quantum vacuum state of each parti- 
cle with oppositely directed momenta and short wavelengths is the adiabatic ground 
state. This is related to the assumption that the universe was in a maximally- 
symmetric state at some time in the past which docs not contain a curvature sin- 
gularity [3]. Due to the accelerated expansion of the universe during infiation, 
quantum fiuctuations are amplified into macroscopical or classical perturbations. 
The early-time fiuctuations lead to the formation of large-scale structure in the 
universe, and also contribute to the temperature anisotropies in the CMBR. For 
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modes whose wavelength is much greater than the Hubble radius, the final state 
for the perturbations is a highly-correlated two-mode squeezed state, with pairs of 
field quanta produced at late times with opposite momenta [4] . Tensor fiuctuations 
in the metric, for example, are predicted to give rise to relic gravitational waves. 
Electromagnetic waves cannot be squeezed during cosmological expansion in the 
same way the tensor (metric) perturbations are, because they do not interact with 
the external gravitational field in the same way as gravitational waves do. 

For the non-rotating, chargeless black hole we shall be considering, the only 
parameter characterising the black hole is its initial mass Mj. In the adiabatic 
approximation for evaporating black holes, for much of the evaporation, frequen- 
cies typically exceed the inverse time- variation scale for the black-hole mass. 
That is, the wave period is much smaller than the timescale of variations of the 
background gravitational field. The black hole interacts negligibly with the emitted 
particles, and the time between successive emissions is comparable with the black 
hole mass [5]. Prior to the final disappearance of the hole, however, particle frequen- 
cies are of the order m, and the frequency of variations of the background space-time 
is comparable to particle frequencies. The particles then interact strongly with the 
evanescent black hole. 

In the cosmological case, a natural length scale is the time-dependent inverse 
Hubble parameter = | is the Friedmann- Robertson- Walker scale factor), 
with the adiabatic approximation k ^ H, or such that the wavelengths of the 
perturbations are less than the Hubble radius aH'^. When the wavelength is com- 
parable with or longer than the Hubble radius, the amplification of the zero-point 
fluctuations takes place. In addition, both the ADM mass and Hubble parameter 
control the redshifting of the radiation in the background space-time. 

Due to the Schrodinger evolution, cosmological perturbations (rotational, density 
and gravitational) in an initial vacuum state are transformed into a highly-squeezed 
vacuum state with many particles having a large variance in their amplitude (par- 
ticle number) and small (squeezed) phase variations. At small wavelengths, the 
squeezing of cosmological perturbations may be suppressed, while it should be 
present at long wavelengths, particularly for gravitational waves [38]. These per- 
turbations also induce the large-angular- scale anisotropics observed in the CMBR. 
Their wavelengths today are of the order of or greater than the Hubble radius. 
The amplification of the initial zero-point fiuctuations gives rise to standing waves 
with a fixed phase, rather than travelling waves. The relic perturbations in the 
high-squeezing or WKB limit can be described as a stochastic collection of standing 
waves. However, it has been suggested that the squeezed-state formalism engenders 
no new physics [24]. 

A prominent theme of this paper is the squeezed-state formalism applied to black- 
hole evaporation, with reference to the comparable inflationary-cosmology scenario. 
Grishchuk and Sidorov [4] flrst achieved a squeezed-state representation for Hawk- 
ing radiation. Their approach, however, referred to final quantum states describing 
the particles escaping to infinity and those falling through the event horizon. In this 
paper, the squeezed-state approach is related to the radiation which remains after 
a black hole has evaporated completely. In an adiabatic approximation, the fixed 
phases correspond to discrete frequencies in the remnant Hawking radiation from 
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the evanescent black hole. This high-squeezing feature, we argue, can potentially 
be observed. The standing waves of the highly-squeezed final state of the Hawking 
radiation originates from choosing to set Dirichlet boundary data on initial and final 
space-like hyp ersurf aces. Such a problem is not well posed, and so we propose com- 
plexifying the proper time interval at spatial infinity in order to obtain a well-posed 
problem. This, we believe, has the effect of avoiding the problems associated with 
the final naked singularity which may be present prior to the final disappearance 
of the black hole. Thus, we conclude that a highly-squeezed final state, which is a 
pure state, is related to the avoidance of the naked singularity, and its associated 
unpredictability. 

After a brief summary of our two-surface method in Section 2, in Section 3 we 
write the quantum amplitude for linearised spin-0, 1, 2 fields in a unified form. Co- 
herent states and generalised coherent states are related to the quantum amplitude 
in Section 4. Section 5 introduces the squeezed-state formalism into the theory via 
a complexification technique. In Section 6, we discuss issues of entropy and classical 
behaviour. Throughout we employ Planckian units G = c = h = 1. 

2 Two-Surface Formulation 

This section is a summary of results which will be elaborated on in ref . [6] . 

Our approach to black-hole evaporation is through a two-surface boundary- value 
problem. We consider the quantum-mechanical decay of a chargeless, non-rotating 
black hole into almost flat space-time and purely outgoing radiation. Data for spin- 
0, |, 1, |, 2 perturbations is set on an initial space-like hypersurface E/, just prior 
to the collapse of an initially diffuse (spherical) 'star' to a black hole, and a final 
space-like hypersurface S^?, sufficiently far to the future of S/ that it catches the 
relic Hawking radiation. The black hole evaporates completely into predominantly 
massless particles, with the final total (ADM) energy equal to the initial mass Mj. 

The classical action S^^''\ i.e., the action evaluated at the solution of the field 
equations, reduces to solely a boundary term, where the boundary comprises S/, 
Ei? and a time- like boundary E°° ({r = Roo}) located at large radius connecting 
E/ to Ei?, the proper time separation at spatial infinity being denoted by T. For 
simplicity, we considered initial data of very weak perturbations, i.e., such that 
the matter and gravitational fields were initially spherically symmetric. The final 
data comprises non-zero perturbations corresponding to the relic radiation after the 
black hole has disappeared. 

Let us specialise to the integer-spin perturbations; the fermions will be considered 
in a future paper. For the wave-like field equations satisfied by the perturbations, the 
'Dirichlet' boundary-value problem on a purely Lorentzian background space-time is 
not well-posed [7]. That is, there may be no solution to the field equations with this 
boundary data, or an infinite number of solutions. Adopting a 'Euclidean' strategy, 
we rotate T into the lower complex plane. Consequently, one might then expect 
there to be a unique classical solution to the nearly-Lorentzian field equations. This 
is just Feynman's ±ie prescription in quantum field theory. 

In a neighbourhood of E/ and E^?, one may employ an adiabatic approxima- 
tion, whereby most of the radiation frequencies exceed the rate of change of the 
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background space-time. One may then 'Fourier' decompose the perturbations into 
standing waves, so that the boundary conditions on E/ and T^p arc satisfied. For 
example, the real linearised massless scalar field A^^\x) which vanishes on S/ and 
is non-zero on is written as 

= If; ^ r dkaoMmAi+{t:r)^^^YU^). (1) 

The background line-element is spherically-symmetric on a time average and has 
the time-dependent Schwarzschild-like form 

where dQ"^ is the two-sphere line-clement. In the limit of wavelengths small com- 
pared to the background radius of curvature, the background geometry at late times 
is approximately a Vaidya space-time. The {^^^(Q)} are spherical harmonics with 
angular momentum / and azimuthal angular momentum m. The {^oki+{t,r)} are 
'radial' functions which are functions of r and slowly-varying functions of t in an 
adiabatic approximation. The {aQkim+} are a set of time-independent amplitudes 
(see below) . Only one set of amplitudes is associated with the final state of the ra- 
diation as we suppose that there is no event horizon since the black hole disappears 
completely. A decomposition similar to eqn.(l) for the spin-1 and spin-2 perturba- 
tions is also possible - see ref. [6]. The boundary- value problem for the fermionic 
fields differs from the bosonic case due to the first-order nature of their actions. 

By taking the limit i?oo oo, one can normalise the radial functions on T^f, 
assuming regularity at the origin {r = 0}, which implies that the radial functions 
are real, and a linear combination of plane waves at infinity. The classical action can 
then be written as a single integral over frequency. For the spin-0, 1, 2 perturbations, 
the classical action has an infinite number of discrete simple poles along the real 
frequency axis at the frequencies a;„ = ^,n = l,2,... This is a manifestation of the 
ill-posedness of the two-surface boundary- value problem. Through our 'Euclidean' 
strategy, we can avoid the poles along the real frequency axis and subsequently 
obtain a probabilistic interpretation for the final boundary data. This is because the 
classical action acquires real and imaginary parts in the complexification procedure. 
At the end of our calculations, wc take the limit that T approaches the real axis, 
a procedure intimately connected with the squeezed-state formalism. In the next 
section, we will compute the quantum amplitude for the integer-spin fields as a 
stepping stone to computing the probability density. 

Additional papers [6] matched our two-surface approach with the method of Bo- 
goliubov coefficients tailored to describing the radiation incoming from and outgoing 
to the null hypersurfaces T~ and X+ respectively. 

3 The Quantum Amplitude 

For the massless spin s = 0, 1, 2 perturbations, we find that the quantum ampli- 
tude or wave functional ^[{osfeimp}; 2^] is given by [6] 

n{askimp}; T] = ATexp {iS^''\{aMmph T]) , (3) 
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where {agkimp} denotes a set of 'Fouricr'-like amplitudes for perturbations with 
spin s, frequency k, angular momentum /, azimuthal angular momentum m and 
polarisation P, and S^'^''^ is the linearised total classical action. (For a massless 
scalar field, the polarisation can refer to the real and imaginary parts of a single 
complex massless scalar field.) Further, is a T-dependent prefactor given below. 
We emphasise that eqn.(3) is just a semi-classical approximation to the full quantum 
amplitude, which is given by a path integral over field configurations which match 
the specified boundary data. For the simplest initial Dirichlet boundary condition 
that the initial perturbations are very weak, we found that 

S^^\{asMmP};T] = --EE E E^^TTX^^/ ^dre'^^^-'^CsimpdtCimpk. 

s l=sm=-lP=± "'O 



- -MiT 
2 ^ 



E (-l)''^^ ., / dkk\zskip\'^\askimP + i-iyPas-kimp\^cot{kT) 



slmP 



(l + s) 



- \miT. (4) 

The final term in eqn.(4) comes from the time-like boundary E°° and gives a plane- 
wave representation for black holes of fixed ADM mass M/ in the absence of per- 
turbations. In eqn.(4), = 2tt,wi = j,W2 = |, and askimP = P{-^y"'a*g _^i _„^p. 
The functions {C,simp(t, satisfy in an adiabatic approximation k ^ ||d — 6| [6] 

e^(f-«)a^(el(''-»)a,e./^p) - d^CsimP + VMt, r)CsimP = 0, (5) 
where Vsip{t,r) is a spin-dependent potential: 



Vsip{t,r) = 



l{l + l) + {l-s') 



2. 2m(t,r) 



(6) 



for s = 0, 1 and odd-parity spin-2 perturbations, and a more complicated term for 
the even-parity metric perturbations which will not be needed in this paper. The 
effect of the back-reaction on the metric in the adiabatic approximation is just to 
replace the Schwarzschild constant mass with the time- and radially- dependent mass 
m{t,r) in Vgip where e^"*^*''') = 1 — ^iMiill. The spin-0 and spin-2 problems arc very 
similar. In cosmology, primordial density perturbations and the relic gravitational 
waves are also described by similar equations. The coefficients {zgkip} relate to the 
boundary conditions at spatial infinity and the regularity conditions at the origin 
for the radial part of the real functions {C,simp{t, r)} [6]. In the first term of eqn.(4), 
one must understand that the limit Roo — cxd is implied. 

Equation (3) can be interpreted as giving a 'coordinate'-representation ampli- 
tude for a set of final field configurations {agkimp} labelled by 'quantum' num- 
bers {sklmP}, given that on the hypersurface E/ {t = 0} the perturbations van- 
ished, i.e., the metric and background matter were spherically symmetric. Hence 
I^Ijc^sfcZmp}; 2^] P is a (conditional) probability density of finding the field in a set 
of final configurations {askimp} a-t asymptotic proper time T. 
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To proceed further, let us discretisc the frequency integral in eqn.(4), with {kj > 
0} the eigenfrequencies of the final radiation located in a spatial volume on — the 
spatial momenta. An explicit expression for the {kj} will not be needed, although 
the continuum limit is recovered in the limit R^o — oo. Denoting j = {slmP}, we 
find that eqn.(3) can be written as (see also Appendix A) 

M{A,};T] ^ iVe-i^^^n 2^3i4^.r) exp (^A%%|^,f cot(%r)) 

OO 

j n=0 

where 

|A,f = 2{-iyws^^^_\z,\'\a, + {-iyPa,,.,^imp\\ (8) 

and En — kj{n+^) is the quantum energy of the linear harmonic oscillator. Note the 
dependence of the quantum amplitude on \Aj\; that is, it is 'spherically symmetric'. 
The {!/„} are Laguerre polynomials with 

4--)(a;)=f:f ^ (9) 

the associated Laguerre polynomials, and L!^\x) — [8]. The {L„(x)} satisfy 

the completeness relation 

oo 

e-"^^Lnix)e--^yLM = d{x - y). (10) 

n=0 

The Laguerre polynomials {L^d^;^)}, where z = x + iy, cannot be written as a 
product of two decoupled wave functions of x and y in an excited state (due to pair 
correlations [22] in a quantum interpretation), but we can write [9] 

Lnix' + y') = E ( ^ ) H,,{x)H,n-2Ay)- (11) 

The wave functional \E'[{Aj}; T] is proportional to a product of delta functions when 
kjT = pn {p being an integer), following from the completeness relation eqn.(lO) - 
see Appendix A. At these focal points, the density of paths - effectively the prefactor 
in the first line of eqn.(7) - diverges. 



4 Coherent States 

The Schrodinger picture wave functions 

*„,(x,-, T) = -e-^^^e-2^^"^L„(x,), (12) 

where Xj = kjAkjlAj]"^, have a strong connection with the exact solution of the 
forced-harmonic-oscillator problem [10]. In this theory, one considers a 1-d harmonic 
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oscillator with mass ji and frequency a;, which is acted on by an external force F{t). 
In this case, the Hamiltonian has the form 

H^\{p' + q') + qF{t). (13) 

Suppose that for to < t < T this force is non- vanishing so that asymptotic states are 
free oscillator states. We want to calculate the probability amplitude Akm to make 
a transition from the free oscillator state \m> (with m particles) at time t' < to, 
before the force begins to act, to the free oscillator state \k> ai time t > T, after 
the force has ceased. Set 

(5 = r dtF{t)e-''^\ (15) 

effectively the 'Fourier' transform of the force. It has been shown that [11] - [13] 

(m > k) 

where A is a real phase and Ak„i is symmetric in k and m. An adiabatic approx- 
imation {z -C 1) indicates that a state which begins as \m > must end up in the 
same state at late times after the time-dependent force has been removed. Then 

Akk = e'^e-"^'Lk{z). (17) 

The probability that there be no change in the number of particles is therefore 
l^fefcP = ^~^{Lk{z)Y ■ Apart from the introduction of mode labels {j} and a neces- 
sary reinterpretation for z, these amplitudes are effectively the semi-classical wave 
functions Eqn.(12) derived from our boundary- value problem. 

A brief derivation of eqn.(16) in the context of the coherent-state representa- 
tion [13] will be useful for the calculations in Section 4.1. Coherent states jo; > can 
be regarded as displaced vacuum states [14]: 

la >= D(q;)|0 >, (18) 

where D{a) — e""^~"*" is a unitary displacement operator {D'^{a) — D~^{a) — 
D{—a)) and the states jo; > are eigenstates of the annihilation operator a with 
complex eigenvalue a. They arc the closest states to classical states in that they 
attain the minimum demanded by the uncertainty principle. Thus, coherent states 
may also be useful when one has a large number of particles, as in the classical 
limit, and when one has some phase information about the state. Coherent states 
form an over-complete set, but are not orthogonal. In terms of Fock number states 
|n>=(^|0>[13]. 



la >=e-2l"l ^ \n > . (19) 
n=o V n\ 
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The coherent state corresponding to a = is the unique state of the oscillator, 
namely, the Fock state \n > with n = 0. Thus, if the system started in a vacuum 
state, the amplitude to then find it in a coherent state |q; > is 



< 0|q; >=< 0|D(q;)|0 >= e-^l°l', (20) 
up to a phase. 

Prom the properties of displacement operators [14] , then 

D{0\a> = DiOD{a)\0> 

= + a)\0 > e^^^'^*-^*"^ 

= le + a > e^(«"*-^*"). (21) 
For later reference, one can show that 

D\a)D(ii)D{a) = L>(At)e">-"''*. (22) 
Using eqns.(19) and (21), 

< m\D{^)\a >= -^(e + a)"^e-^(H^+l€P+2ra). (23) 
Also from eqn.(19), 

oo n 

< m\D{0\a >= e-^l"l' V ^ < m\D{C)\n > . (24) 

n=o V n! 

Equating eqns.(23) and (24) then gives 

oo / 1 

(1 + y)"*e-^l«l' = e^l^l' V n^r"""?/" < ^1^(01^ > ■ (25) 



n=0 » 



The generating function for the associated Laguerre polynomials [8] 

oo 

(1 + yre-y^ = y: \y\ < i, (26) 



n=0 

then implies that the matrix element between initial and final states can be written 
as 

ml 



< m\D{0\n >= p. r-"e-tl^l^Lf-)(|er), (27) 



which is eqn.(16) up to an unimportant phase factor. Interchanging m and n gives 
the same result. 

4.1 Generalised Coherent States 

There is an interpretation for these amplitudes in terms of generalised coherent 
states \n,a > of the harmonic oscillator [12]. Defining 

\n, a >= e~'^"^D{a(t))\n >, (28) 
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then in the Fock representation 



\n,a >= ^ < m\D{a{Q))\n > \m > e '^""K (29) 



m=0 



For generahsed coherent states (see eqn.(33) below), the ground state {n — 0) is 
a coherent state and not a vacuum state. Generahsed coherent states are to the 
coherent states what the Fock states \n > are to the vacuum state, that is, excited 
coherent states. In addition, one has 

I — — [ (Pa\n,a >< n,a\, (30) 

TT J 

<n,p\n,a> = L„(|a - /^He^'^-^d^l'+l/^l'), (31) 
<n,P\^> = — / (i^Q;L„(|Q;-/3|^)e^*"e-^l"l' < n,Q;|V' >, (32) 

TT J 

for an arbitrary state \ip >, where 

— J (fa = J d{Rea)d{Ima) , 
and / is the identity operator. From eqn.(31) with P — 0, then 

<n,0\n,a> = <n\n,a> ^ e"^l"I^L„(|Q;|^), (33) 

again giving eqn.(12) up to a phase and normahsation. The initial state may be 
seen not as a vacuum state, but as a Fock state, and the final state as a generahsed 
coherent state with the same n. 

An interesting interpretation for the amplitudes eqn.(27) is that these are the 
matrix elements for the transition from the state |A; > to |m > under the influence of 
a gravitational wave [15] , with the force F{t) proportional to the Riemann curvature 
tensor component Rxtxt{t)'- 

F{t) = ficHR,t.t{t) = —ficHd^hlJ, (34) 

where I is the distance between two particles each of mass |/x along the x-axis, h^J 
is the transverse tracclcss gravitational wave component [16] and x is the change in 
separation of the masses. In the context of black-hole evaporation, one can under- 
stand the time-dependent force as being active during the time-dependent phase of 
the collapse to a black hole and its subsequent complete evaporation. The duration 
of the 'force', therefore, is comparable to the lifetime of the black hole. A space- 
time would begin in an (almost) static state in the far past, where the perturbation 
modes were effectively free, pass through an intermediate time-dependent phase, to 
end up in a static or, rather, quasi-stationary state in the far future, where again 
the perturbation modes are free. In this 'sandwich'-space-time picture, particles are 
created by the space-time curvature, or non-adiabatic behaviour of the fields in a 
time-dependent metric. What the above calculations indicate is an explicit math- 
ematical connection between the theory of forced harmonic oscillators and certain 
amplitudes relating to the dynamical evolution of black holes. 
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5 Squeezed-State Formalism 

In this section, we see how by rotating the proper time at infinity T into the complex 
plane, and with spherically symmetric initial matter and gravitational fields, one 
obtains a quantum-mechanical squeezed-state interpretation for the final state of 
the Hawking radiation. 

Grishchuk and Sidorov [4] were the first to explicitly formulate particle creation 
in strong gravitational fields in terms of squeezed-states, although the formalism 
does appear in Parker's original paper on cosmological particle production [18]. In 
ref. [4], it was shown that relic gravitons (as well as other perturbations), created 
from zero-point quantum fluctuations as the universe evolves, should now be in 
a strongly squeezed state. Squeezing is just the quantum process of parametric 
amplification. In ref. [24] however, the claim was made that there was no new 
physics in employing the squeezed-state formalism. 

5.1 Squeezed States - Introduction 

A general one-mode squeezed state (or squeezed coherent state) is defined as [17] 

\a,z>^ Si{r,ip)D{a)\Q >= Si{z)D{a)\Q >= Si{z)\a >, (35) 

where D{a) is the single-mode displacement operator, 5*1 (r, ip) is the unitary squeez- 
ing operator {S\{z)Si{z) = Si{z)Sl{z) = 1) for \a,z >: 

S,{r,ip) = Si{z) = e-^^^^'-'*''''\ (36) 

where z = re^^*'^. Further, Si{r, Lp)Si{r' , ip) — Si{r + r' ,(fi). With a\a >— a\a >, 
then setting b = SiaSl implies that b\a,z >= a\a,z > and b = acoshr -|- 
a^e^**^ sinhr. The state |q;,2 > is a Gaussian packet displaced from the origin in 
position and momentum space. While the real squeezing parameter r > deter- 
mines the magnitude of the squeezing, the squeezing angle </'(— f<(/7<f), gives 
the distribution of the squeezing between conjugate variables. The operators a and 
are annihilation and creation operators, respectively. The squeezed vacuum state 
occurs when a = 0: 

\z>=\0,z>^Si{z)\0> . (37) 

The high-squeezing limit corresponds to r ^ 1, where the state \z > is highly 
localised in momentum space. The state with r = is the ground state. 

Consider the amplitude A =< a, z\D{fi)\a, z >=< z\D''{a)D{fi)D{a)\z >. One 
canuseeqn.(22) to show that A =< z\D{fi)\z > e^*-^'"^"*^) =< 0\D{r])\0 > e2»-f™("», 
where rj = /^.coshr — /i*e^"^ sinhr, and then use /i = |yu|e*'^ to show that [15] 

|^|2 _ g-|r;P _ g-|Mp[cosh2r-cos2((/)+(p)sinh2r] 

which is independent of the initial displacement a. Equation (38) will appear when 
we match our space-like hypersurface approach with the method of Bogoliubov 

coefficients [6]. 

Single-mode squeezed operators do not conserve momentum since they describe 
the creation of particle pairs with momentum k. Two-mode squeezed operators, 
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(38) 



however, describe the creation and annihilation of two particles (waves) with equal 
and opposite momenta. A two-mode squeeze operator has the form [29] 

S^{t, ^) = e^(e--a+— e--aVl)^ ^ ^-1 ^ g^iv, ^ + |), (39) 

where a±, a± are annihilation and creation operators for the two modes, respectively. 
One can also introduce a unitary rotation operator 

R{e) = e-^^(»+"+ +«-«-), (40) 

with < ^ < 27r, R{9)R{9') = R{9 + 9'), R\0 >= |0 >,RaR^ = e'^a. However, 
rotation does not influence particle creation [29] . 

The variances of conjugate operators p and q are given by Aq = q— < q >, Ap = 
p— < p > . In the squeezing formahsm, Aq and Ap differ greatly, while they are 
equal and the minimum possible for coherent states, deemed to be the most classical 
of quantum states. The name 'squeezed' refers to the fact that the variance of one 
variable in a conjugate pair can go below the minimum allowed by the uncertainty 
principle (the squeezed variable), while the variance of the conjugate variable can 
exceed the minimum value allowed (the superfluctuant variable) [22] [26] [27]. The 
superfluctuant variable is amplified by the squeezing process, and so it is possible to 
observe macroscopically, while the subfluctuant variable is squeezed and becomes 
unobscrvable. In particle production, for example, the number operator is a su- 
perfluctuant variable, while the phase is squeezed. Squeezed states are essentially 
purely quantum-mechanical in origin. 

Hawking radiation in the squeezed-state representation was first discussed in 
ref. [4]. The squeeze parameter vj was related to frequency and the black-hole 
mass through tanhr^ = e~^'"'^^^ ,9j = ipj (see also refs. [21] and [32]). The vacuum 
quantum state in a black-hole space-time for each mode is a two-mode squeezed 
vacuum. 

5.2 Analytic Continuation and the Leirge-Squeezing Limit 

Consider the Schrodinger picture quantum state eqn.(7) associated with the final 
hypersurface Ep, and define 

mA^y.T] = J{2is\n{k,Tm{A,];T] = iVe-^^^^n^xp V-Ak,k,\A,\' coi{k,T) 

j 3 

= N exY> {^S^''\{A,};T]) . (41) 

In order for our two-surface boundary-value problem to be well posed, we argue that 
T must be a complex quantity. This is because for real T, the 'sum' in eqn.(41) 
diverges due to the simple poles on the real-frequency axis at /c ^ = cr„ = = 
1,2,..., assuming that /cjjAjp remains finite and non-zero near kj — an- At the 
frequencies there may be an infinite number of classical paths (solutions) 

joining the initial and final data, or perhaps none at all. Should a solution exist, it 
may not depend smoothly on the boundary values. These are features of ill-posed 
boundary- value problems. Yet, if T is deformed into the complex plane, the poles 
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along the real-kj axis would be displaced into the complex-/c plane. The complex 
boundary- value problem itself guarantees that the classical solution is complex ana- 
lytic (strong ellipticity) . One then has a reasonable expectation that there will now 
be a smooth, unique classical solution of the field equations joining the initial and 
final data. Essentially, this is the ie prescription of quantum-field theory. In gen- 
eral, therefore, our background metric 7^j, (and hence g^„) is not a real Lorentzian 
metric but perhaps a strongly-elliptic metric. This may permit the path integral 
over the boundary data to be well-defined. Even for a complexified T, one still has 
a time, \T\, whose 'gradient' is everywhere future-directed and time-like. By con- 
sidering the dimensionless quantity T/2Mj, one also sees that analytic continuation 
of the complex quantity T is equivalent to adding a small positive imaginary part 
to the mass, as is the case in the Damour-Ruffini tunnelling approach to black-hole 
evaporation [19]. 

In terms of the late stages of black-hole evolution, we may use analytic continu- 
ation to tackle the crucial question of whether future null infinity X+ is a Cauchy 
hypersurface for space-time. If we model the evaporation process with a classical 
Lorentzian space-time metric, it is well-known that a momentarily-naked singular- 
ity is inevitable [20]. Cosmic censorship, a postulate of classical general relativity, 
is seemingly transcended through the quantum-mechanical evaporation of the black 
hole. Consequently, a Bogoliubov transformation between initial and final modes 
would be forbidden. There would then be no one-to-one correspondence between 
pre-coUapse configurations and post-evaporation configurations, no S-matrix and, 
therefore, there will be evolution from a pure initial state to a final density matrix. 
Our reasoning is that a non-singular asymptotically-fiat initial state on a suitable 
Cauchy hypersurface can never evolve to a final (even momentarily) naked singu- 
larity if we complexify the proper-time separation at spatial infinity, so that the 
future is entirely predictable from the past. If measurements at I'^ are insufficient 
to determine the state at T~ , one needs boundary conditions ~ or another set of 
'coordinates'- on the Cauchy horizon which forms after the black hole disappears. 
Boundary conditions along the Cauchy horizon represent constraints on the initial 
conditions as configurations on a final Cauchy hypersurface must be sufficient to 
determine the state on an initial Cauchy hypersurface. That there is no singularity 
or edge to space-time is a necessity if the path integral for quantum gravity is to be 
a no-boundary wave function [3]. 

In the presence of a time-independent Schwarzschild black hole, say, one com- 
monly 'Wick' rotates the real time t 90° into the lower half- plane {t — —it) [1]. 
One now has a Euclidean time separation at spatial infinity and the space-time 
boundary is real and positive-definite. Where there is a future event horizon, one 
has the additional boundary component at the axis r = 2Mq. In the Euclidean 
regime, a conical singularity is present at r = 2Mo unless the Euclidean time is 
periodically identified with period SttMo [2]. Thus, the coordinate singularity at 
r = 2Mo as well as the curvature singularity at r = is avoided, since r > 2Mo for 
Euclidean signature. The manifold for which r > 2Mq and < r < SttMo is the real 
Euclidean section of the Schwarzschild solution. Boundary conditions for physical 
fields then need not be specified at the curvature singularity, which does not he on 
the Euchdean section, and these fields are also regular at the horizon. This was 



14 



Euclidean time 



Lorentzian time 




Figure 1: Infinitesimal deformation of asymptotic proper- time separation by an angle S into the 
lower complex plane. 



interpreted as quantum cosmic censorship, and with the regularity at r = 2Mo, one 
is seemingly summing over all possible configurations which lie inside the black hole 
r < 2Mq. 

Hence, let us write (see fig.l) 



T = \T\e-'^, 0<d < 



TT 



(42) 



Below we shall consider the case of infinitesimal S. Then, from eqns.(41) and (42) 

mA,};T] = $[{A,};|T|,5] 



_/Vg- iiM^ |r| cos 5g- iMj |r| sin i J-j- 



--Akjkj\Aj\'^ coth{kj\T\ sin 5 - iipj{\T\,6)) 



(43) 



where 
and 



nRj{\T\,S)^ 



sinh(2A;j T 


sin 6) 


2[cosh^(A;j|T 


sin 5) — cos^ ipj] 



(Pj{\T\,S) = —kj\T\ cos S, 

, nij{\T\,s) 



(44) 



sin 2ipj 



2[cosh {kj\T\ sin 5) — cos^ ipj] 



(45) 



15 



One can also write eqn.(43) as 

mAj}; \T\,S] = iVe-i^^^l^l<=°^'5e-^^^l^l^''^^nexp 



2 ■'Ml- e^^^nanhr, ' ' ^ 



where we have formally set 



tanhrj(|T|,5) = e 



-2kj\T\sinS 




Therefore, 

= tanh(A;^|T| sin 6). (48) 

We recognise eqn.(46) as the coordinate-space representation of a quantum-mechanical 
squeezed state [30] [32], with rj{\T\,S) the squeeze parameter and ipj{\T\,5) the 
squeeze angle. The evolution of the squeezed state is taken into account by the 
|T| dependence in Vj and (fj, which are in general complicated functions of time. 
Equation (44) is more familiar in the limit of infinitesimal S, whence neglecting 
0((5^) terms, Lpj{\T\,5 <^ |) ~ —kj\T\, corresponding to free evolution. The high- 



frequency limit kj\T\ ^ 1 corresponds to 



for sin S 0. 



Computing the probability density |$[{74j}; |r|, 5] p, we find for small 5 



m{Aj};\n6]\'^\N 



2^-Mi\T\5 



Y[exp 



COth Ej 



1 + 



^^^^^Akjkj\Aj\'^ 



where < Ej = kj\T\S < 1, 
Then, from eqn.(48) 



i.e., < 5 < {kjfy 



with kj 



2Mikj, T 



Ej^e 



-2r,- 



Sj < 1, 



(49) 



2Mj ■ 

(50) 



corresponding to rj ^ 1, which is the high- squeezing limit. 

In the squeezed- state formalism, this limit is deemed to be the classical limit when 
the average number of particles in the final state is large: < Nj >= sinh^ rj ~ ^e^^^ 
for ^ 1. Here, the limit of infinitesimal 5 is also a quantum-to-classical transition, 
from a Euclidean (5 ~ |) theory, where there is no concept of time, to a semi- 
classical Lorentzian theory, with infinitesimal d, where the notion of a classical time 
parameter appears. Such issues are related to the quantum cosmology program of 
Hartle and Hawking [3] in the emergence of the classical universe from a smooth 
Euclidean (quantum) origin, where there is no cosmological-singularity. 

Indeed, the semi-classicality of the state eqn.(41) is intimately related to the 
high-squeezing limit [24]. If we just consider Mj as a fixed parameter in the theory 
and not a functional of the final field configurations {xj}, then the WKB condition 
is met when 







sin 2(/9j 






sinh(2A;j T sin 5) 



sin 2(fj sinh 2rj | ^ 1, 



(51) 



which is satisfied in the high-squeezing limit Ej -C 1 even if sin2<^j = 0. The fi- 
nal state of the remnant Hawking flux, therefore, becomes more classical in the 
WKB sense in the limit S ^ 0. In this limit, one can effectively consider the 
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final perturbations as being represented by a classical probability distribution func- 
tion [4] [24] [40]. As in the inflationary scenario, the perturbations on the black-hole 
background space-time, which had a quantum-mechanical origin, cannot be distin- 
guished from classical stochastic perturbations, without the need of an environment 
for decoherence. It is the field amphtudes which one observes, and these amphtudes 
have a classical phase-space distribution. The squeezed nature of the final quantum 
state for the stochastic gravitational-wave background produced by inflation is not 
stationary [4]. 

The initial conditions for the perturbations in the black-hole case also have an 
inflationary analogue. In cosmology, the assumption is that at some early 'time' 
just prior to inflation, the modes are in their adiabatic ground state. This originates 
from the assumption that the universe was in a maximally-symmetric state at some 
time in the past [3]. A similar assumption obtains in our black-hole case, where 
we assumed that the initial perturbations were very weak so that the initial matter 
and its gravitational fleld were spherically symmetric. 

We assume that T ^ 1; that is, we observe the black hole at inflnity at times 
much greater than its collapse timescale, which is of order ttM/ [16]. If kj is mod- 
erately large, then for small Sj (keeping 0{ej) terms) 

m{Aj};\T\,S]\' 

where a;„ = Acj„ = -j^ and we used kj > with — > as kj ^ to 
convert the sum Yl,^=-oa i'^to Hence, assuming the interchange of j and n 

sums in the continuum limit for the {kj} frequencies is valid, then 

oo 

|$[{yl,}; |T|,5 ^ 0]|2 = |iV|2 [] [] g-Aa;„a,„|A,„,^Pp^ ^53) 

slmP n=l 

SO the frequencies {a;„} dominate the flnal state. This is the same as in ref. [6] 
where we used contour integration to obtain the probability density. Equation 
(52) describes a Gaussian non-stationary process in that the variance is an oscilla- 
tory function of the asymptotic proper time. Rather than travelling waves, one is 
now dealing with standing bosonic waves, where the amplitudes for left- and right- 
moving waves are large and almost equal, similar to the cosmological scenario [4] - 
see eqn.(l). These standing waves imply a correlation between particles with oppo- 
site frequencies (and azimuthal angular momentum m) in the flnal state, indicative 
of a pure state. We suggest that the discrete frequencies {a;„} are, in principle, 
observable feature of the relic Hawking radiation, either directly or indirectly. 

The presence of the delta function in eqn.(52) indicates that in the high-squeezing 
limit, the random variable ipj associated with the flnal state is squeezed to discrete 
values independently of the quantum numbers {slmP}. Note that it is only the 
squeeze phases {^j} of the perturbations which are flxed and correlated in the 
high-squeezing limit. There are phases in our theory, however, which we deem to 



|^|2g-M,|T|5-Qg^p 



"4 + sin^(fc,|T|)^^^'^^''^^''' 



[TVl^JJexp 



—/S.kjkj\Aj'\^ ^ /S.UJn5{kj — UJn) 



n=l 



(52) 
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have a random and even distribution; namely, the phases associated with the ini- 
tial state of the fluctuations incoming from T~, and the phases associated with the 
spatial field distribution labelled by {agkimp}- The final field amplitudes for the 
remnant Hawking radiation reaching the null hypersurface X"*" {bj} can be matched 
on Ei? with the amplitudes {aj} associated with the space-like hypersurface [6]. 
The {bj} are in turn connected, through a Bogoliubov transformation, with ampli- 
tudes {cj'} for travelling waves incoming from X^. A positive-definite probability 
distribution can be obtained by averaging over the random phase of the {cj/} fluc- 
tuations. Therefore, it is only the phases {(pj} associated with the standing waves 
which do not have a random character. 

The normalisation factor | A^p is determined by integrating over the dimensionless 
variables {yj}, where xj = AkjkjlAjl"^ = Mfyj, so that the sum of all probabilities 
of all possible conflgurations {yj} is unity. Treating the initial ADM mass Mj as a 
fixed parameter rather than as a functional of {xj}, we find that 

Mf coth Ej 



w = n 



^ _^ sm2(fc,|T|) 



sinh'^ 



= JJ[cosh 2rj - cos 2ipj sinh 2rj]"^Mf (54) 
j 

In fact, in the high-squeezing hmit 5 — > 0, final states which differ by a coordinate- 
dependent phase e~2^^''^ are deemed to be physically equivalent. If we considered 
the ADM mass as a functional of {xj}, then the normalisation factor would be 
altered for finite 6, as would the value of the entropy - see below. Thus, there is an 
ambiguity in the normalisation factor. In addition, this is related to the ambiguity 
in the form of the entropy in the high-squeezing limit with and without a surface 
term included in the Lagrangian [30]. In the limit Sj 0, then 



M{A,};\T\,S^Qf 



slmP n=l 



This equation is reminiscent of the Wigner function corresponding to the vacuum 
state in the limit of infinite squeezing [40], giving the classical 'trajectories' of the 
system. The Wigner function cannot in general be interpreted as a classical proba- 
bility density for finite squeezing, that is, finite 6, except for Gaussian states. 

A consequence of the high-squeezing behaviour is that the variance in the am- 
plitudes {xj} is large, so that there are large statistical deviations of the observable 
power spectrum from its expected value. This is just a manifestation of the uncer- 
tainty principle. Indeed, with respect to the first expression in eqn.(52), and using 
eqn.(54), we find that 



< yj >i*|2 £j 



n=l 

= [PjiMj)]-\ (56) 

where <>[$[2 denotes the expectation value with respect to the probability density 
It is reasonable to suggest that pj{Mj) represents a microcanonical density-of- 
states. 
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In inflationary cosmology, the oscillation phases of standing waves have fixed 
values, giving rise to zeros in the power spectrum, which are characteristic of the 
CMBR. In this case, the power spectrum of the cosmological perturbations in the 
present universe is not a smooth function of frequency. The standing-wave pattern, 
due to squeezing, induces oscillations in the power spectrum. This in turn produces 
Sakharov oscillations [24] [42] , produced by metric and scalar perturbations, in the 
distribution of higher-order multipoles /(/ + 1)C/ of the angular correlation function 
for the temperature anisotropics [38] [39] in the CMBR for all perturbations at a 
fixed time whose wavelength is of the order of or greater than the Hubble radius 
defined at this time. That is, the peaks and troughs of the angular power spectrum 
have a close relationship with the maxima and zeros of the metric power spectrum. 
However, for long wavelengths, the power spectrum is sufficiently smooth. Sakharov 
oscillations, therefore, exist as a result of the squeezed nature of the scalar and 
metric perturbations. 

6 Entropy and Squeezing 

There have been many accounts of how to determine the entropy generation in 
the squeezing formalism [22] - [30]. Hu and Pavon [31] were the first to associate 
entropy generation with the monotonic increase in the average particle number with 
time, induced by parametric amplification in a vacuum cosmological space-time. As 
squeezing is the quantum analogue of parametric amplification, one would expect 
that the squeezed-state formalism can compute entropy production. This is indeed 
the case, although, as with any entropy calculation, the nature of the coarse-graining 
must be specified. For squeezing, this is particularly relevant as squeezed evolution 
is unitary, i.e., there is no loss of information in principle in the evolution of the 
initial pure state to the final pure squeezed quantum state. How one chooses to 
measure the observables associated with the final squeezed state determines the 
entropy. One can reduce the final density matrix with respect to a Fock or coherent 
state basis [30], or use eigenstates of the superfluctuant variable [22] [26] [27]. In 
refs. [22] [26] [27], the loss of information comes from the increased dispersion of 
the superfluctuant operator. 

In the classical limit of large average particle number, corresponding to the large- 
squeezing regime ^ 1, a universal form for the entropy density growth ASj 
obtains for each mode [22] - [30]: 

AS J ~ 2rj, rj » 1, (57) 

irrespective of the particular coarse-graining. In particular, on averaging over the 
squeezing angle (pj of each Fourier mode, one obtains eqn.(57) [22] [23] [25] - [28]. 

One can determine whether coarse-graining with respect to the squeeze angle is 
appropriate from the prescription discussed in ref. [32]. Calculating the entropy S 
from eqn.(49) using eqn.(54), then 

S = -llldyjP(yj)\nP(yj) 
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= 1 + ^ In [M72(e2^^ sin^ + e-^^^' cos^ ^j)] . (58) 
j 

Thus, the entropy eqn. (58) comes from our ignorance of the final radiation config- 
uration. The constant term is not important. In the high-squeezing limit 

5~ 1 - ^lnM| + 2^r^ + ^lnsinVj-. (59) 

3 3 3 

Evidently, sm^pj — for particular values of the frequency. Even if sm. (pj — > 0, 
for Tj oo and ignoring constant terms, then one may argue that ^Sj ~ 2rj. 
However, an alternative derivation keeping 0(£|) terms gives from eqn. (54) 

oo 

3 "=1 

= l-^lnp,(M,). (60) 

3 

How we take the high-squeezing limit is therefore important. The effect of the 
final term in eqn. (59) is to reduce the entropy from the maximal value given in 
eqn. (57). In the cosmological scenario, the primordial gravitational- wave back- 
ground entropy is significantly smaller than eqn. (57) due to the information about 
the initial inflationary state of the universe manifesting itself in the primordial peaks 
in the multipole spectra of the CMB temperature anisotropy [43]. The presence of 
these peaks, therefore, is incompatible with a totally random squeezing angle, and 
so a coarse-graining with respect to pj is not possible. A similar conclusion may 
be argued in our case. We will comment on this in another paper [6], where we 
discuss the discreteness of the radiation frequencies in the context of Bekenstein's 
discrete event-horizon-area theory, and how information about the initial state may 
be contained in the spectral lines. 



6.1 Classical Predictions 

We now discuss how strong peaks in the wave function lead to definite predictions. 
In quantum cosmology, wave functions are commonly peaked about correlations 
between coordinates and momenta. In rcf. [34], identification of such correlations 
came via the Wigner function. If the Wigncr function W{p^ q) factorised into a 
function of position q and a function of momentum p, then no correlation between 
q and p is predicted. When W{p, q) is peaked about some region in phase space 
p = f{q), then the wave function predicts this particular correlation. 

An alternative proposal for measuring correlations was given in ref. [35]. Here, 
projection onto coherent states, where momentum and position are equally known, 
as in classical theory, is employed for predicting classical correlations from a general 
Wigner function. It was shown that in the harmonic-oscillator case, the correlation 
between p and q was such that the Hamiltonian equalled the classical energy. We 
arrive at similar conclusions. In ref. [34], the predicted correlation was between 
the Hamiltonian equated to the quantum energy. The problem with this approach 
was that crude approximate Wigner functions bore no resemblance to the exact 
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Wigner functions. How Wigncr functions appear in the context of our two-surface 
formalism will be considered in another paper. 

If the quantum state ^' is 'sufficiently' peaked about a region in the phase 
space, we observe correlations between the observables which characterise this re- 
gion. Phase-space configurations for which ^ is small are precluded and will not 
be observed. Where ^ is neither small nor sufficiently peaked, no predictions can 
be made [36] . These conclusions are drawn from a new interpretation of quantum 
mechanics where wave functions are not associated with probabilities. 

With the ADM mass Mj as a parameter, we look for predictions from the Heisen- 
berg picture wave functional "^l^^Xj) — ^e~^^i Ln^ixj). Were we to restore our 
units, then the argument of the Laguerre polynomial would be proportional to TT^. 
Thus, in the limit of large argument, we use for large x [8] 

Lk{x) t^. (61) 

Prom a dimensional argument, Xj oc —f ^1. In this case also, therefore, the 

pi 

approximation in eqn.(61) applies. One can now find the peak in the wave function 
as a function of Xj at 

nj^lkjAkj\Aj\'n-\ (62) 
Considering the spin-0 case, for example, from eqn.(8) and eqn.(62) 

Uj = 2'!rkjAkj\zj\'^\aj + {-lyPag^^kjimpl'^- (63) 
In a previous paper [6] , we showed that for spin-0 perturbations 

= 27ik,\zj\^\aj + {-lYPa,,_k,imp\^ (64) 

where the {bj} are Fourier amplitudes associated with the radiation reaching X"*"; i.e., 
eqn. (64) matches the positive- frequency decomposition for massless spin-0 particles 
reaching X+ (travelling waves) with our two-surface method for field configurations 
on the {t — T} hypersurface Hp (standing waves). Hence, rij — Akj\bj\'^. Thus we 
find that 

Y,hn,k,=Y.Ak,k,\b,\'^Mi. (65) 
j j 

The right-hand side is just the final total energy in the massless spin-0 fiuctua- 
tions, which equals the initial ADM mass Mj, as demanded by our boundary- value 
problem. The left-hand side is just the total energy of the radiated particles. 

7 Conclusion 

In this paper, we have illustrated many aspects of the two-surface formulation for 
linearised integer-spin fields propagating in an evaporating black-hole space-time. 
When the proper-time separation T at infinity between the initial and final hy- 
persurfaces is deformed infinitesimally into the lower complex plane, and one has 
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spherically-symmetric initial fields, one obtains a quantum-mechanical squcczcd- 
state formalism. The large squeezing limit is equivalent to the WKB limit and 
corresponds to an infinitesimal angle in the lower complex-T plane. This we believe 
is related to the emergence of time in a semi-classical Lorentzian space-time from 
a timeless Euclidean regime. As the highly-squeezed final state is a pure state, our 
complexification technique has seemingly avoided the unpredictability associated 
with the momentarily-naked singularity, which is believed to be present prior to the 
complete disappearance of the black hole. 

As in the cosmological scenario, we found that the bosonic perturbations on 
the black-hole background can be deemed to be a stochastic collection of standing 
waves rather than travelling waves in the high-squeezing limit. This leads to the 
prediction of peaks in the power spectrum of the relic Hawking radiation analogous 
to the Sakharov oscillations in the CMBR. 

Many of the features discussed in this paper may be valid in the instance of non- 
spherical initial fields; i.e., non-zero initial perturbations. In this case, there will 
be cross terms in the classical action involving the initial and final field configura- 
tions. One can subsequently diagonalise the classical action resulting in a product 
of squeezed states, and so the results in this paper are valid. For rotating black 
holes, similar conclusions to those arrived at in this paper arc possible. One must 
also include the fermionic perturbations, which are described by anti-commuting 
Grassmann variables. 
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Appendix A — Derivation of the wave function eqn.(7) 

To arrive at eqn.(7), we set a — e"*'^^^, y — and x — Xj in Mehler's formula [9] 




p=0 



(66) 



where 




(67) 



are Hermite polynomials which satisfy 




(68) 



Then, 



oo 




(69) 



p=0 
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where Ep — kj{p + ^) and 



U^^) - , (70) 



I.e. 

|2 



da;j|^p(a;j)| =1. 
Only even terms contribute to the sum in eqn.(69) as [8] 

^2p+i(0) = 0. (72) 

Hence, 

oo 

[27risin(^,T)]-^e^^^?^°*('=^^) = Y.^-'^"'^ ^•.^{x,)^,^^- (73) 
In addition, set Xj — (^)^i?e(Aj), = (^)^7m(Aj). Then, 

^^ = E e-^(^^^+^v)^^,^(^,)^,^,(^,)^,^(0)^,,,(0) 

oo p 

= V'2p'(a;j)V'2p-2p'(%)V'2p'(0)V'2p-2p'(0) 
p=0 p'=0 

where y"-*^ = Akj. Inthehmit |T| — > 0+ (or rather kj\T\ <^ 1), as the eigenfunctions 
ipp{xj) form a complete orthonormal set 

lim [27risin(%T)]-^e*^^l^il'-t('=.^) = 5(a;,)5(7/,) 

^ <5(^H(^)^^.)- (74) 

This suggests that Aj — > as kj\T\ 0+, agreeing with our initial conditions of 
very weak perturbations. We now use the identity eqn.(13). Taking the product 
over all j and introducing a normalisation factor A^, then 



j p=0 \^ J 

which is eqn.(7), where we have included the contribution from the time-like bound- 
ary. One can confirm that eqn.(75) gives us the first part of eqn.(7) by simply using 
the generating function for Laguerre polynomials [8]. 
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